Numerical simulation on compressible turbulence by spectral method

The numerical and physical issues of simulations on compressible turbulence are reviewed in the present paper. An outline of the global spectral methods and the progress of recent local spectral methods are illustrated. Several typical subjects in this field are studied, including homogeneous isotropic turbulence, autoignition in premixed turbulence, interaction between flames and turbulence, and shock wave in turbulence. The results of the numerical simulations are discussed, enabling us to discover and to understand the physical phenomena which have not been solved by experiments.     

高黏性流动有限元模拟的迭代稳定分步算法

分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于 某个临界值方能使算法稳定. 然而在高黏性流动模拟中，已有的显式和半隐式分步算法由于 其显式本质，必须采用小时间步长计算，不但降低了计算效率，同时也常与为使分步算法稳 分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于 某个临界值方能使算法稳定. 然而在高黏性流动模拟中，已有的显式和半隐式分步算法由于 其显式本质，必须采用小时间步长计算，不但降低了计算效率，同时也常与为使分步算法稳 定必须满足的最小时间步长要求冲突. 本文目的是构造一种含迭代格式的分步算法，它能在 保证精度的前提下大幅度地增大时间步长. 方腔流和平面Poisseuille流数值计算结果证实 了此特点，该方法被有效应用于充填流动过程的数值模拟.     

Numerical analysis of pressure oscillations over axisymmetric spiked blunt bodies at Mach 6.80

The pressure oscillations over a forward facing spike attached to an axisymmetric blunt body are simulated by solving time-dependent compressible Navier–Stokes equations. The governing fluid flow equations are discretized in spatial coordinates employing a finite volume approach which reduces the equations to semidiscretized ordinary differential equations. Temporal integration is performed using the two-stage Runge–Kutta time stepping scheme. A global time step is used to obtain a time-accurate numerical solution. The numerical computation is carried out for a freestream Mach number of 6.80 and for spike length to hemispherical diameter ratios of 0.5, 1.0 and 2.0. The flow features around the spiked blunt body are characterized by a conical shock wave emanating from the spike tip, a region of separated flow in front of the hemispherical cap, and the resulting reattachment shock wave. Comparisons of the numerical results are made with the available experimental results, such as schlieren pictures and the surface pressure distribution along the spiked blunt body. They are found to be in good agreement. Spectral analysis of the computed pressure oscillations are performed employing fast Fourier transforms. The surface pressure oscillations over the spike and phase plots exhibit a behaviour analogous to that of the Van der Pol equation for a self-sustained oscillatory flow. Received 28 February 2001 / Accepted 17 January 2002     

Direct simulation of breakdown to turbulence following oblique instability waves in a supersonic boundary layer

The late stages of transition to turbulence in a Mach two boundary layer are investigated by direct numerical simulation of the compressible Navier-Stokes equations. The primary instability at this Mach number consists of oblique waves, which are known to form a pattern of quasi-streamwise vortices. It is found that breakdown does not follow immediately from these vortices, which decay in intensity. The generation of new vortices is observed by following the evolution of the pressure and vorticity in the simulation, and analysed by consideration of vorticity stretching. It is found that the slight inclined and skewed nature of the quasi-streamwise vortices leads to a production of oppositely signed streamwise vorticity, which serves as a strong localised forcing of the shear layer alongside the original vortices, formed by convection and stretching of spanwise vorticity. The shear layer rolls up into many new vortices, and is followed by a sharp increase in the energy of higher frequencies and in the skin friction.     

可压缩材料平面应变滑移线场的数值解和边界条件

简要介绍了可压缩材料平面应变滑移线理论。导出了滑移线场应力和速度数值解。处理了应力边界条件，阐明了应力和速度间断规律。举例说明了边界条件处理和数值解步骤。     

Swirl effects on variable-density jet mixing in confined flows

An experimental investigation on swirl effects on inhomogeneous confined jet mixing in a combustor configuration is reported. The confined swirling flow was simulated by a swirler with a central jet mounted in a cyclindrical tube. Helium and air jets set at different velocities were injected into the confined swirling air flow. The resulting flow fields due to two vane swirlers with constant vane angles of 35° and 66° were compared. Results show that the 35° vane swirler produces a solid-body rotation core with a slope about twice that created by the 66° vane swirler. It is the behavior of this solid-body rotation core that determines jet mixing rather than the swirler vane angle. Consequently, the coaxial jet decays much faster, the mixing is more intense, and the turbulence intensities are higher for the 35° vane swirler. In view of these results, combustor designers should be more concerned with behavior of the solid-body rotation core produced by the swirler, instead of the swirler vane angle.     

三维N-S方程计算隐式有限体积法

利用NND有限差分格式，发展了一种新的完全隐式的有限体积数值方法，以求解与时间相关的N－S方程．对通过单元体界面的无粘流和粘性流通量均作隐式处理．对绕流钝锥体和不同攻角的气动辅助实验飞行器的高超声速粘性流和化学反应流获得了定常数值解．对流加热率和流场电子密度的计算值与实验数据进行了比较，符合较好，证实了本方法的精确性．     

快速求解叶栅流动的边界元方法

本文利用边界元方法，通过ｋｒｉｃｈｈｏｆｆ变换将描述叶栅流动的控制方程转换成线性方程。并将广义ｋ－Ｊ条件与边界积分方程联立求解，避免了非线性项和叶片出气角的迭代计算。完成了一种快速求解任意迥转面叶栅流场的计算程序，实用表明与其它数值方法及实验结果符合较好，具有快带、简明、实用的特点。     

Numerical modeling of the disturbances of the separated flow in a rounded compression corner

Numerical modeling of the time-dependent supersonic flow over a compression corner with different roundness radii is performed on the basis of the solution of the two-dimensional Navier-Stokes equations in the regimes corresponding to local boundary layer separation. The development of unstable disturbances generated by local periodic injection/suction in the preseparated boundary layer is calculated. The results are compared with those of similar calculations for a flat plate. It is shown that the natural oscillations of the boundary-layer second mode stabilize in the separation zone and grow intensely downstream of the reattachment point. The acoustic modes excited within a separation bubble are studied using numerical calculations and an asymptotic analysis.     

Analysis and Modeling of the Turbulent Diffusion of Turbulent Kinetic Energy in Natural Convection

Buoyant flows often contain regions with unstable and stable thermal stratification from which counter gradient turbulent fluxes are resulting, e.g. fluxes of heat or of any turbulence quantity. Basing on investigations in meteorology an improvement in the standard gradient-diffusion model for turbulent diffusion of turbulent kinetic energy is discussed. The two closure terms of the turbulent diffusion, the velocity-fluctuation triple correlation and the velocity-pressure fluctuation correlation, are investigated based on Direct Numerical Simulation (DNS) data for an internally heated fluid layer and for Rayleigh–Bénard convection. As a result it is decided to extend the standard gradient-diffusion model for the turbulent energy diffusion by modeling its closure terms separately. Coupling of two models leads to an extended RANS model for the turbulent energy diffusion. The involved closure term, the turbulent diffusion of heat flux, is studied based on its transport equation. This results in a buoyancy-extended version of the Daly and Harlow model. The models for all closure terms and for the turbulent energy diffusion are validated with the help of DNS data for internally heated fluid layers with Prandtl number Pr = 7 and for Rayleigh–Bénard convection with Pr = 0.71. It is found that the buoyancy-extended diffusion model which involves also a transport equation for the variance of the vertical velocity fluctuation gives improved turbulent energy diffusion data for the combined case with local stable and unstable stratification and that it allows for the required counter gradient energy flux.